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ABSTRACT 

Achieving a robust determination of the gas density profile in cluster outskirts is a crucial 
k>" ■ point in order to measure their baryonic content and to use them as cosmological probes. 

The difficulty in obtaining this measurement lies not only in the low surface brightness of the 
ICM, but also in the inhomogeneities of the gas associated to clumps, asymmetries and ac- 
cretion patterns. Using a set of hydrodynamical simulations of 62 galaxy clusters and groups 
we study this kind of inhomogeneities, focusing on the ones on the large scale that, unlike 
C*~) . clumps, are the most difficult to identify. To this purpose we introduce the concept of residual 

dumpiness, C-jz, that quantifies the large-scale inhomogeneity of the ICM. After showing that 
this quantity can be robustly defined for relaxed systems, we characterize how it varies with 
radius, mass and dynamical state of the halo. Most importantly, we observe that it introduces 
an overestimate in the determination of the density profile from the X-ray emission, which 
translates into a systematic overestimate of 6 (12) per cent in the measurement of Af gas at 
i?20o for our relaxed (perturbed) cluster sample. At the same time, the increase of C-r with 
\ radius introduces also a ^2 per cent systematic underestimate in the measurement of the 

hydrostatic-equilibrium mass (Mh e ), which adds to the previous one generating a systematic 
~8.5 per cent overestimate in / gas in our relaxed sample. Since the residual dumpiness of 
the ICM is not directly observable, we study its correlation with the azimuthal scatter in the 
X-ray surface brightness of the halo, a quantity that is well -constrained by current measure- 
ments, and in the ^/-parameter profiles that is at reach of the forthcoming SZ experiments. 
We find that their correlation is highly significant (r$ = 0.6 — 0.7), allowing to define the 
azimuthal scatter measured in the X-ray surface brightness profile and in the y-parameter as 
robust proxies of C-n- After providing a function that connects the two quantities, we obtain 
that correcting the observed gas density profiles using the azimuthal scatter eliminates the bias 
in the measurement of M gas for relaxed objects, which becomes ± 2 per cent up to 2i?200> 
and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the 
systematics on the measurements of Af n c and / gas , although a significant halo to halo scatter 
remains. 



Key words: Cosmology: large-scale structure of Universe - X-rays: galaxies: clusters - meth- 
ods: N-body simulations. 
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1 INTRODUCTION 

Clusters of galaxies are the most massive virialized systems of the 
Universe, they form in the knots of the cosmic web from which 
they continuously accrete material in the form of dark matter, gas 
and galaxies. Their importance for astrophysics is crucial, since 
they enclose information o n the large-scale structure formation (see 
lKravtsov~& Borgani 2 0121 . for a review) and provide constraints 
on cosmological parameters which are independent with respect 
to cosmic microwave back ground (CMB), type-la supernovae and 
galaxy surveys (see, e.g.. I Allen et al.l201ll and references therein). 

The precision with which clusters can be used as cosmolog- 
ical probes depends on the accuracy in the measurements of their 
mass, not only in terms of the total mass but also for its baryonic 
fraction /b = Mb/M to t- If the value of /b measured inside clus- 
ters matches the cosmic one (flb/fi m ) it can be used to obtain an 
independent estimate of f2 m . However we know that the assump- 
tion of /b — fib /f2 m holds only at large radii thus requiring the 
achievement of accurate observations also of cluster outskirts. In 
the recent years this point raised a great interest of the astrophysi- 
cal community towards the study of clusters around and beyond the 
virial radius. 

When moving outside the core the intracluster medium (ICM) 
is often affected by the presence of inhomogeneities and sub- 
structures, whose impact on the observed properties has been 
studied with a statistical approach (see, e. g., Jeltem a"et alj|2005l ; 
iBQhringer etaill20ld lAndrade-Santos et al.ll2012l) . The difficulty 
in observing cluster outskirts lies in the very low X-ray surface 
brightness of the ICM, that drops below the diffuse extra galactic X- 
ray background at r ^i^orFksee lRoncarelli et al.l2006al lbl). A pos- 
sible m ethod consists in sta cking the images of different objects, as 
done bv lEckert et alj d2012l) who measured the density profiles and 
the azimuthal scatter in the X-ray surface brightness of 3 1 ROS AT- 
PSPC objects, showing a clear segregation between Cool-Core 
(CC) and No-Cool-Core (NCC) systems over the radial range 0. 1- 
0.8 Rwo, with the first ones exhibiting a scatter of 20-30 per cent, 
which corresponds to density variations of the order of 10 per cent. 
Above 0.8 -R200, the azimuthal scatter increases up to values of 60- 
80 per cent both in CC and NCC clusters, suggesting that they expe- 
rience common physical conditions in shaping their X-ray surface 
brightness profiles. However the limit of the stacking technique is 
that it does not allow a precise determination of the temperature 
profiles of single objects. In this framework a great step forward 
has come from the Suzaku X-ray satellite that, thanks to its very low 
instrumental background, provided the first robust spectroscopical 
analys es of the regions c l ose to -B200 for the brightest objects (see , 
e.g.. iGeorge et al l 120091 ; iBautz et al.1 120091: iHoshino et alj |2010|: 
Kawaharad a et alj|20ld : ISimionescu et alj l201ll : iHumphrev et al.1 
2012l : ISato et alj|2012r) . These works also highlighted the presence 



of a substantial azimuthal scatter, likely associated to clumps or 
inhomogeneities on the large scales, induced by major or even 
minor mergers th at can cause sloshing and swirling motions in 
the ICM (see alsolSanders & FabiarJl2012l; Ichurazov etai1l2012l : 



ISimionescu et alj|2012h 

A promising complementary view with respect to X- 
rays can come from t he observations of the thermal Sunyaev- 
Zeldovich effect (SZ , Sunvaev & ZeldovicrJ Il970l : iBirkinshawl 
1 19991 : iRephaeli et alj|2005h . The SZ signal is associated to tem- 
perature fluctuations in the CMB spectrum that are directly pro- 



1 In this work -R200 is defined as the radius enclosing an average density 
equal to 200 times the critical density of the Universe p c = . 



portional to the integral of the electron pressure along the line of 
sight. The SZ effect, in combination with the X- ray emission, can 
therefore probe the triaxial structure of the ga s (|Pe Filippis et al.1 
120051: iMorandi etai]|2012l ; ISereno et ai]|2012l) . Moreover, if the 
source is sufficiently extended to be resolved spatially even with 
the large beams available with present instruments such as Planck 
(as for nearby, hot systems like Coma, IPlanck Collaborationll2012l) 
or multi-pixel b olometer arrays (as the APEX-SZ experiment, 
iBasu et al]|2010l or Boloc am at the Caltech Submillimeter Obser- 
vatory, Savers et"al1l2012h it can be used also to directly infer the 
pressure profile, under the assumption of the spherical symmetry of 
the ICM distribution, out to a significant fractio n of the virial radius 
(see also lWalker et al.ll2012tlEckert et al.ll2013l) . 

From the theoretical point of view the modelisation of clus- 
ter outskirts is affected by different uncertainties with respect to 
cluster cores. On one side different feedback mechanisms do not 
affect significantly the behavior of temperature and density pro- 
files, with gravity that cons titutes the dominating physical pro- 
cess teoncarelli etai] [2006b). On the other hand, the presence of 
shocks and turbule nce may lead to the break of the hydrostatic equi- 
librium (see, e.g. , Iapichino & Niemeverll2008l; IVazza et al.|[2009l : 
iBurns et"al]|2010l : iNagai & Laull201 ll . and references therein). In 
addition to that, hydrodynamical simulation s showed how the outer 
ICM is affected by t he presence of clumps dRoncarelli et al.l2 006b; 
Nagai & La ull201rl) and inhomogeneities dKawahara et alj 2008 ; 
Vazza et al.1 1201 If) that may bias the reconstruction of the clus- 



ters' properties and explain the observed a zimuthal scatter (see also 
IVazza etalj|2012l ;l Zhur avleva et al These large-scale inho- 

mogeneities, both in gas densities and temperatures, are responsible 
for half of the ~ 30 per cent underestimate in the X-ray recon- 
structed hy drostatic mass, a s highlighted by numerical sim ulations 
(see, e.g jRasia etal]|2012l;|piffaretti & Valdarninill2008l) . The re- 
maining part is mostly due to the residual bulk motions of the ICM, 
that preserve energy in non-thermalised form, not traceable with the 
measures of the density and temperature profiles alone ( see, e.g., 
iRasia et al]2006r . lBattagiia et alJl2012UNelson et alj2012l and ref- 
erences therein), but potentially detectable with spatially-resolved 
X-ray spectroscopy (s ee, e.g., the works on the Coma clusters by 
ISchuecker etai1l2004l and lChurazov et alj|2012l) . 

In the work presented here we study the gas inhomogeneities 
present in the outskirts (close to ifeoo and beyond) of galaxy clus- 
ters and groups using a set of hydrodynamical simulations, with the 
objective of analyzing how the different physical processes may af- 
fect the degree of density fluctuations, together with the halo mass 
and its dynamical status. We concentrate our efforts in the defini- 
tion of the density inhomogeneities associated to large-scale fluctu- 
ations, for which we introduce the concept of residual dumpiness, 
i. e. the dumpiness of the ICM after obvious clumpy structures 
have been removed. We discuss how this phenomenon may bias 
high the X-ray density profile measurements and explore, for the 
first time, the possibility of obtaining a direct measurement of this 
quantity via the analysis of the observed azimuthal scatter in the 
X-ray and {/-parameter profiles of single objects. This method of 
estimating the residual dumpiness of the ICM allows us to propose 
a technique to reduce this bias and, therefore, to improve the pre- 
cision of the measurement of the total mass and of the gas mass of 
galaxy clusters. 

This paper is organized as follows. In the next Section we de- 
scribe the set of hydrodynamical simulations used for our study. 
In Section [3] we characterize the gas inhomogeneities of our sim- 
ulated haloes, we provide the definition of the residual dumpiness 
(C-r.) and we study its dependence on cluster properties. Section [4] 
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presents our results on the correlation between the azimuthal scat- 
ter and large-scale inhomogeneities of cluster outskirts. In Section[5] 
we quantify how the residual dumpiness biases the measurements 
of the gas mass, the hydrostatic-equilibrium mass and the gas frac- 
tion of the haloes and present a method to correct them using the 
observed azimuthal scatter. We summarize and draw our conclu- 
sions in Section[6] 



2 THE SIMULATED CLUSTERS 

The clusters object of our work belong to a set of 29 high-resolution 
re-simulations of galaxy cluster regions. A detaile d description of 
the wh ole procedure to obtain them can be found in lBonafede et al.l 
d201ll) . Here we provide a summary. 

2.1 Simulation parameters 

The cosmology assumed is a flat ACDM model with f2 m = 0.24 
and Qa si — f2 m = 0.76, a Hubble parameter h — 0.72 (be- 
ing the Hubble constant Ho = 100 ft km s _1 Mpc -1 ); we fix 
the primordial power spectrum of the DM fluctuations with slope 
n — 0.96 and normalization coherent with erg = 0.8. 

With these assumptions we carried out a large dark-matter 
(DM) only simulati on using the G ADGET- 3 code, an evolved ver- 
sion of GADGET- 2 dSpringelll2005l) . with a periodic box of 1 h~ x 
Gpc on a side and then identified the clusters in the z = out- 
put with a Friends of Friends (FOF) algorithm, with linking length 
fixed to 0.17 times the mean inter-particle separation. The La- 
grangian regions around 7i? v ir from the centre of the 24 most 
massive objects (all with M v - lv > 10 15 /i~ Mq) were identi- 
fied, traced back to their initial positions and then re-simulated 
at higher resolution using the zoomed initial conditions technique 
(Torm erTet al.l 1997h . These re- simulations were run adding the hy- 
drodynamical part by turning on the Smoothed Particles Hydro- 
dynamics (SPH) code implemented in GADGET-2. This was done 
assuming fib = 0.045, with mass resolutions of 8.4 x 10 8 Ii~ 1 Mq 
and 1.6 x 10 8 /i _1 Mq for the DM and gas particles, respectively, 
and with a Plummer-equivalent softening length for the gravita- 
tional force fixed to £ = 5ft" 1 kpc in physical units at z < 2 
and kept fixed in comoving units at earlier times. 

In order to add statistics at the lower mass scales, another set 
of 5 galaxy clusters with masses M200 = 1 — 5 x 10 14 h~ 1 Mo was 
selected at z — and re-simulated with the same method, reaching 
a total of 29 re-simulated regions. 

2.2 Feedback models 

The set of 29 re-simulations was run with three different physi- 
cal implementations. More detailed descript ions of these physical 
models can be found in lPlanelles et alj a2012h . 

(i) NR: non-radiative runs. They do not include any physical pro- 
cess except gravitation and hydrodynamics. Although unrealistic, 
they are useful as a test to check the impact of the various physical 
mechanisms. 

( ii) CSF: runs including cooling, star formation, metal enrichment 
and galactic winds. The radiative coolin g rates are computed fol- 
lowing the procedure of Iwiersma et al. I (l2009h . considering also 
the contribution of metals in the hypothesis of a gas in (photo- 
)ionis ation equilibrium by using the CLOUDY code jFerland et alj 
1998). Star formation is followed according to the multi-phase 



model of Springel & Hernquist (2003]) which also includes en- 
ergy release by superaovae (SNe). Other feedback sources include 
heating from a sp atially uniform, time dependent UV background 
(Haardt&Madau] |l99d) and metals released by t ype-II and type- 
la SN e and AGB stars according to the model of iTornatore et al.l 
J2007I) . Finally, a kinetic feedback associated to galactic ejecta is 
implemented assuming a mass upload proportional to the star for- 
mation rate and a wind speed of « w = 500 km s~ . 
(Hi) AGN: like CSF, but with wind speed reduced to v w = 350 
km s _1 and including also the feedback associated to gas accretion 
onto supermassiv e black h o les (B Hs). Its implementation is sim- 
ilar to the one of ISpringell (120051) . with the inflow of matter that 
proceeds at Bondi rate up to the Eddington limit. In our implemen- 
tation SPH particles stochastically selected for providing accretion 
into the BH feed it only with 1/4 of their mass at time, thus mim- 
icking a more continuous flow of material. The amount of kinetic 
energy released by each BH particle is given by 

E{ eed = e r £fM| H , (1) 

being e r and et the radiative efficiency and the fraction of energy 
coupled to the gas, respectively. In our implementation these two 
free parameters were fixed at e r = 0.1 and et = 0.05, increasing 
to ef = 0.2 when the accretion r ate is smaller than one-hundredth 
of the Eddington limit (see also Isiiacki et alj|2007l ; iFabian et alj 
2010). We will take this model as our reference one. 

2.3 Definition of the halo sample 

Since the volume of each of the 29 re-simulations encloses a spher- 
ical region of radius much larger than the virial radius of the main 
object, it is common to find other haloes included in the final 
2 = outputs. Hence, in order to increase the statistical signif- 
icance of our results and to characterize better also the low mass 
objects, we included in our sample all the 'secondary' haloes with 
A/200 > 10 14 ft _1 M0 in at least one of the three runs. The values 
of A/200 were computed starting from the FOF catalogues and ap- 
plying a Spherical Overdensity algorithm. This increases our halo 
sample to a total of 62 objects in different mass ranges. A histogram 
of their masses is shown in Fig.[TJ 

In order to study the dependence of the dumpiness on the 
mass of the halo and on its dynamical properties, we divide our 
objects into subsamples: given the amplitude of our initial sample, 
this can be done without losing statistical robustness. In particular 
we will adopt the following two definitions: 

- clusters! groups (28 and 34 objects, respectively), according to 
their M200 being higher/lower than 5 x 10 14 /i _1 M , 

- relaxed! perturbed (30 and 32 objects, respectively), accord- 
ing to a definition based on the dumpiness and described in Sec- 
tionl3~3l 

The histogram of Fig. Q] shows also these classifications. It is 
worth to mention that clusters tend to be slightly more perturbed 
(16 and 1 2 relaxed) than groups (16 and 18) due to their more active 
dynamical state. 



3 DENSITY INHOMOGENEITIES 
3.1 The gas dumpiness 

A commonly used indicator of the degree of inhomogeneity of a 
medium is its dumpiness. Even if it is normally used to measure 
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Figure 1. Distribution of M200 for the sample of 62 simulated haloes in our reference (AGN) run. Cluster names correspond to the ones given in Bonafed e et alj 



( 20 11), with the addition of D25-D29 that indicate the low mass objects re-simulations. The main haloes of each re-simulation are dubbed DXX-a (bold face) 
while their satellites are indicated with letters from -b to -e in decreasing mass order. The vertical dashed line indicates the mass limit of 5 X 10 14 h~ 1 Mq 
used to separate groups (34 objects) and clusters (28). Objects marked in red correspond to perturbed (32) haloes, while black indicates relaxed (30) ones, as 
defined in Section|3~3"l 



the amount of small clumps in an approximately uniform medium, 
its general definition allows us to use it also to quantify the amount 
of inhomogeneity associated to large-scale density fluctuations. 

The dumpiness or clumping factor C of a fluid element is 
defined by the following formula: 



C 



< p >v 



(2) 



< P >v 

where p is the fluid density and the brackets <>v indicate the 
average calculated over its volume. The dumpiness is therefore de- 
fined as being equal to unity for a perfectly uniform medium and 
> 1 otherwise, with higher values indicating higher levels of inho- 
mogeneities. When writing explicitly the integrals, eq. ((2} can be 
expressed as as 

V 



C : 



M 2 



p z dV 



(3) 



where M and V are the mass and volume of the gas element, re- 
spectively. Since we are using SPH simulations, for a given distance 
r from the cluster centre, we compute C(r) by converting eq. ((3) 
into 



Or) 



(4) 



where the sum Y2i is computed over the particles with distance 



between r and r + dr, im and pi are the mass and density, respec- 
tively, of the i-th SPH particle and V shc ii = f 7r[(r + dr) 3 - r 3 ] is 
the volume of the shell used to compute the quantity. 

For the purpose of our work, for every halo we compute the 
dumpiness and the other physical quantities in the range < r < 
2R200, in 50 equally spaced bins. 



3.2 Separating small clumps from large-scale 
inhomogeneities: the residual dumpiness 

The amount of dumpiness of the ICM is caused by two different 
phenomena: the presence of small dense clumps and the density 
fluctuations on larger scales associated to asymmetries in the large- 
scale accretion pattern, with the first one that constitutes the domi- 
nating contribution. Here we describe a method that proves useful 
to isolate these two components in order to treat them separately. 
This aspect is particularly important when we consider that the 
physical properties of clumps and their abundance depend on the 
processes of cooling and star formation and, consequently, their 
presence in our simulated clusters is subject to the uncertainties 
associated to the implementation of physical mode ls in hydrody- 
namical codes. Moreove r, several works (see, e.g., iMitchell et al.1 
120091 ; ISiiacki et alj|2012l and references therein) highlighted how 
SPH simulations produce a higher quantity of dense clumps with 
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respect to Eulerian ones. This is associated to the different numeri- 
cal viscosity, intrinsic in the codes themselves, that causes a lower 
dissipation efficiency in SPH simulations. 

To this purpose, we fo ll ow the volume-selection scheme de- 
scribed in Ron carelli et al.l d2006bl) to compute the profiles of 
galaxy clusters close to and beyond R200. This method consists 
in sorting the SPH particles belonging to a given radial bin accord- 
ing to their physical density. Then, starting from the most diffuse 
one, we sum up their volumes (defined as Vi = mi/pi) until we 
reach a fixed fraction, 99 per cent in our case, of the total volume of 
the radial bin and consider only these particles to compute C: the 
remaining particles are identified as clumps and considered sepa- 
rately in our computation. This procedure proved to be effective in 
excising all the dense and cold clumps that are highly dependent 
on the physical assumptions, as well as the small X-ray bright re- 
gions that are usually masked out in observational analyses, thus 
providing a good description of the global properties of the ICM 
such as dens ity, temperature and X-r ay surface brightness (see the 
discussion in lRoncarelli et alj2006bl) . Recently, a similar approach 
has al so been applied to Eulerian simulations by Zhuravleva et al. 
d2013h who approximate the density distribution of the gas in dif- 
ferent radial bins with a lognormal distribution and mark as clumps 
the cells with density above values fixed in terms of the a of the 
distribution (see also Khed ekaretafll2012h . In fact, although ap- 
plied in completely different numerical simulations (grid-based and 
SPH), the two methods are formally identical, with their threshold 
of fcut = 3.5 that matches our 99 per cent criterion. For what 
concerns the possible biases caused by unresolved clumpy gas, we 
have verified that they have a small impact on our results. For more 
details about this point and about the clumps physical properties, 
we refer to the Appendix [A] We use this method to define the resid- 
ual dumpiness associated to large-scale perturbations, Ctz, and to 
separate it from the total one: this allows us to note that the value of 
C-jz is about one order of magnitude lower with respect to C inside 
R200 (see the discussion in Section |33t . 

We have also checked that in most of the cases lowering the 
threshold value to 95 per cent (i.e. increasing by a factor of 5 the 
volume of gas considered as small-scale clumps) produces negli- 
gible changes in Ctz. However, when the cluster is perturbed by 
the presence of a large halo (e.g. an infalling group), lowering the 
density threshold for the clumps identification can change the final 
result by up to 50 per cent in the radial bins corresponding to the 
halo distance. A sketch of the application of this method is shown in 
Fig.[2]for two different cases. These X-ray surface brightness maps 
show that when cutting the 1 per cent densest volume (from top to 
middle panels) we are able to remove all obvious clumpy structures 
that are associated to the brightest peaks. When cutting the 5 per 
cent of the volume (bottom panels) the surface brightness remains 
almost unchanged for a relaxed halo (left) while for a dynamically 
perturbed one (right) the infalling halo on the bottom right is pro- 
gressively removed, indicating that it is difficult to define a precise 
threshold to separate the two components. The effect on the defini- 
tion of Ctz can be clearly seen in the plot of Fig. [3] where we show 
the residual dumpiness as a function of radius computed for these 
two haloes by adopting the two different volume threshold: while 
for D17-a the value of Ctz is almost independent on the threshold 
chosen up to -R200, for D12-a the presence of the infalling halo at 
r ~ 0.5i?200 produces a difference of about 0.25. 

For this reason we use the difference between the value of Ctz 
computed with these two limits to introduce another useful halo 
classification. We sort our haloes according to the maximum rela- 
tive difference between C 99 (r) and C 95 (r) for r < R200 and split 
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Figure 3. Residual dumpiness as a function of distance from the centre 
computed for the two haloes shown in the maps of Fig. ff] D17-a (black 
lines, lower values) and D12-a (red lines, higher values). The value of 
Ctz has been computed adopting two different thresholds for the volume- 
clipping method: 99 per cent (solid lines) and 95 per cent (dashed lines). 
The physical model assumed is the reference one (AGN). 

our 62 haloes roughly into two equal subsamples: we define a halo 
as relaxed when this difference is less than 8 per cent, and per- 
turbed otherwise. In this way, we end up with 32 perturbed haloes 
and 30 relaxed ones. For this last set we can consider our value of 
Ctz to represent a robust estimate of the amount of inhomogeneities 
associated to large-scale asymmetries. We also verified that our 
classification has a good correspondence with an observational-like 
classification based on X-ray images. We refer the reader to Ap- 
pendix|B]for the details on this point. 

3.3 Dependence on physics and environment 

We show on the top panel of Fig. [4] the results of the dumpiness 
as a function of radius for our whole sample of objects, computed 
considering the three physical implementations described in Sec- 
tion |Z2] When radiative cooling is included (CSF and AGN mod- 
els) the level of dumpiness is very high and ranges from ~ 3 close 
to the centre up to ~ 10 close to 7?200- When extending to outer re- 
gions the value of C increases exponentially reaching values of the 
order of 100 at 2i?2on, with no significant difference when includ- 
ing BH feedback. When neglecting gas cooling (NR), the values 
drop down by a factor of 5-10 over the whole range, indicating that 
the value of C is mainly associated to small cold clumps. 

However, when comparing our results with observational es- 
timates one must consider that the cold dense gas at T < 10 5 
K does not produce any observable dumpiness in X-ray images. 
This makes the plotted NR profile model a reference for the ob- 
served value of C without considering any detection of emitting 
cl umps. For this reason we c onclude that the value of ~16 obtained 
by ISimionescu et alj d201ll) for the Perseus cluster is not realistic, 
since it is sufficient to exclude cold gas from our simulations to 
obtain significantly lower values. 

The lower panel of Fig.|4]shows the values of Ctz, which cor- 
responds to the value of the dumpiness computed after applying 
the 1 per cent volume-clipping described in Section POl With this 
method the difference associated to the physical implementations 
disappears almost completely, even for the NR model: this indi- 
cates that large-scale inhomogeneities do not depend on the phys- 





Figure 2. Bolometric (0.1-10 keV) X-ray flux (in units of counts s _1 cm -2 ) of a relaxed cluster (D17-a, left column) and a perturbed one (D12-a, right 
column). The size of each map is I-R200 P er side. In the first row all the particles have been considered, in the second and third ones the volume cut has been 
applied considering the 99 per cent and 95 per cent threshold, respectively. 
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Figure 4. Upper panel: gas dumpiness of our simulated haloes as a function 
of distance from the centre, for the AGN (solid line), CSF (dashed) and the 
NR model (dot-dashed). The values represent the median of 62 objects, the 
grey-shaded region encloses the quartiles of the AGN model (the quartile 
regions of the other two runs have similar sizes). Lower panel: same as 
upper panel but for the residual dumpiness, i.e. after applying our volume- 
clipping method. 



Figure 5. Upper panel: residual gas dumpiness, i.e. after applying the 
volume-clipping method, for clusters (solid line) and groups (dashed line) 
as a function of distance from the centre for our reference model. The val- 
ues represent the median of the sample (28 and 34 objects, respectively), 
and the grey-shaded region encloses the quartiles of the clusters sample. 
Lower panel: same as upper panel but with relaxed (30 objects, solid line) 
and perturbed (32 objects, dashed line) haloes. 



ical process that occur in the ICM, but on the intrinsic dynamical 
properties of the halo itself. In this case the median of Cm com- 
puted on the whole sample has a value very close to 1 (e.g. almost 
uniform medium) close to the centre and grows constantly up to 
1.3-1.4 at i?2oo- There is a significant variance between the differ- 
ent haloes with about 25 per cent of objects having values of ~ 1.4. 
Outside Rqoo the value of Cm reaches ~ 2 at 2.R200 . together with 
an increase of the dispersion between the objects: this reflects the 
intrinsic difference between the environment of the clusters' and 
groups' outskirts that can contain or not infalling haloes and ac- 
creting filaments. 

This independence with respect to the phys ical implementa- 
tion is partially in contrast with what obtained by Zhuravleva et al. 
J2013h that observe a slightly higher degree of inhomogeneity of 
the bulk in their NR runs, of the order of AC = 0.1 - 0.2. In 
fact, also in our simulations we see a small (AC ~ 0.04) sys- 
tematic excess in our NR model, associated to the hot-dense tail 
of the ICM distribution. This discrepancy probably indicates a dif- 
ferent cooling efficienc y in the simulations analysed here and by 
IZhuravleva et al.l d2013b . 

It is interesting to see how the trend of Cm varies when con- 



sidering the mass of the halo and its dynamic status. We show in 
the top panel of Fig. [5] the results of Cm computed when separat- 
ing the haloes according to their mass into groups and clusters (see 
Section [2~3l . Clusters have on average a higher value of Cm of the 
order of ~ 0.1, with larger differences outside -R200: this is due to 
the fact that massive objects have higher probability of accreting 
material even at later epochs, while galaxy groups are dynamically 
older, more dynamically stable and, therefore, more uniform. When 
separating relaxed and perturbed objects according to our defini- 
tion (see Section [3~2t the difference is even more remarkable, as it 
can be seen in the bottom panel of Fig. [5] While close to the cen- 
tre the two samples have very similar values of Cm, the presence 
of accreting structures pushes the residual dumpiness of perturbed 
objects beyond 1.5 at i?200- On the other side, relaxed haloes show 
all the same trend with Cm raising linearly to 1.3-1.5 at -R200: the 
residual dumpiness of these objects should be considered as an in- 
dicator of the amount of inhomogeneities common to all haloes, 
even the more dynamically stable, and connected to their departure 
from spherical symmetry and to the density fluctuations associated 
to the large-scale accretion patterns. 

We observe as well that the differences in Cm associated to 



8 M. Roncarelli et al. 



the clusters/groups classification can be explained completely by 
the different correlation with the relaxed/perturbed one, (i. e. clus- 
ters have higher values of C-r. because they are on average more 
perturbed than groups), making the latter classification the more 
interesting for our purposes. We note also that the trends th at we 
observe are similar to those found bv lZhuravleva et al1(l2013l) . (we 
refer to their 'bulk CSF' results), both in terms of radial dependence 
and with respect to the dynamical state of the haloes. 



4 THE AZIMUTHAL SCATTER AS A DIAGNOSTIC OF 
THE RESIDUAL CLUMPINESS 

4.1 The azimuthal scatter 

As said, the dumpiness of the ICM is not a directly measurable 
quantity. Here we investigate the possibility of obtaining an esti- 
mate from the observed azimuthal scatter of the X-ray and SZ pro- 
files of the haloes. 

Following Vazz a et al. we define the azimuthal scatter 

of an observable quantity x as 



\ 



N ^ 



5(r) 



(5) 



where N is the number of azimuthal sectors, Xi(r) is the radial 
profile of the quantity in a given sector and x(r) is the average 
over the TV sectors at distance r from the centre. In our case we 
fix TV = 12, so that each sector corresponds to 30 degrees: this 
value is enough to account for all of the main azimuthal fluctua- 
tions a ssociated to typic al ICM inhomogeneities (see the discus- 
sion in lVazza et alj201 lh . For the purpose of our work we consider 
the scatter in the X-ray surface brightness (computed in different 
bands) and in the thermal SZ effect. 

We compute the X-ray surface brightness for every halo in the 
same 50 radial bins and in the 12 different sectors by assuming 
an APEC emission model dSmithetal.ll200lh . by fixing the red- 
of our clusters to z = 0. Since our AGN model follows 
also th e enrichment of metals with the recipe of iTornatore et al.l 
d2007t) . we consider also the contribution from different chemi- 
cal species in the computation of the emissivity of the SPH par- 
ticles (the details of the procedure are described in lRoncarelli et al] 
l2012h . For what concerns the SZ scatter, the method to compute the 
va lue of the ^/-parameter from our SPH simulation is the same as 
in lRoncarelli et al] d2007l) . For these computations we are consid- 
ering all the SPH particles, without applying the volume-clipping 
method (see Section [3~2t , 

We show in the upper panel of Fig. [6] the azimuthal scatter 
as a function of distance from the centre for the value of the y- 
parameter and the surface brightness in three different X-ray bands. 
At each distance we show the median computed over the whole 
sample of 62 haloes in the AGN model. The general observed trend 
is an increase of the azimuthal scatter with radius, with the value 
of cr y going from ~0 at the centre to ~0.3 at ifeoo- The scatter 
associated to the X-ray surface brightness is 2-3 times higher at all 
distances, with lower energy bands (0.3-0.8 and 0.5-2 keV) having 
higher values with respect to the hard 2-10 keV band, indicating 
that these inhomogeneities are associated to gas at temperatures 



(0.3-0.8 keV) 
(0.5-2 keV) 
(2-10 keV) 





Figure 6. Upper panel: azimuthal scatter of the y-parameter (solid line) and 
of the X-ray surface brightness in the 0.3-0.8, 0.5-2 and 2-10 keV bands 
(dashed, dotted and dot-dashed line, respectively) as a function of distance 
from the cluster centre. The values represent the median computed over the 
whole sample of 62 haloes in 12 azimuthal sectors. Lower panel: same as 
upper panel but for the 2-10 keV surface brightness only, computed with 
three different physical implementations: AGN (solid line), CSF (dashed) 
and NR (dot-dashed). The grey-shaded area encloses the quartiles of the 
clusters sample for our reference model. 



around 0.5-1 keV, that has the peak of its emission at soft X-ray 
energies. 

It is interesting to note that when studying how the azimuthal 
scatter varies with respect to the assumed physics and to the halo 
classifications, the trend is very similar with respect to C-r., rather 
than the global clumping factor C (refer to Section \3.3\ . For in- 
stance, the lower panel of Fig. [6] shows that the value of itsb is 
almost independent of the assumed physical model. When looking 
at the dependence on the mass and the dynamical status of the halo 
(upper and lower panels of Fig. [7] respectively) we can see that on 
average galaxy clusters show a slightly higher azimuthal scatter, 
while perturbed systems have values of a which are higher by a 
factor of 2 with respect to relaxed ones at i?200- 



2 Since we are interested in the azimuthal scatter of the X-ray surface 
brightness and not in its absolute value, fixing the redshift has effect only in 
the definition of the X-ray bands. 



4.2 Correlating the dumpiness with the azimuthal scatter 

Since we have shown that dumpiness and azimuthal scatter have 
similar radial trends, here we address directly the question whether 
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Figure 7. Upper panel: azimuthal scatter of the 2-10 keV surface bright- 
ness for clusters (solid line) and groups (dashed) as a function of projected 
distance from the centre for our reference model. The values represent the 
median of the two samples (28 and 34 objects, respectively), and the grey- 
shaded region encloses the quartiles of the clusters sample. Lower panel: 
same as upper panel but with relaxed (30 objects, solid line) and perturbed 
(32 objects, dashed line) haloes. 



a tight correlation exists between these variables and if there exists 
a possible relation that can connect one with the other. 

To this purpose we perform a statistical analysis of these vari- 
ables by computing the Spearman's rank correlation, rs, of C and 
of Ctz with the scatter in the X-ray and SZ profiles of our clusters. 
We restrict our analysis to the interval 0.15 < r/i?20o< 1.5, in 
order to cut out the core, and to the sample of the relaxed haloes for 
which the definition of Ctz is more robust. Given the large sample 
of haloes available, these restrictions do not affect the robustness of 
our computation. The results are shown in Table [TJ 

The main conclusion that can be drawn from these values is 
that the azimuthal scatter of both X-ray and SZ profiles has a very 
high degree of correlation with the values of Ctz, i- e- with the inho- 
mogeneities on the large scales: in fact, the values of rs = 0.6—0.7 
indicate that it is possible to describe Ctz as a monotonic increas- 
ing function of a y and osb- On the other side when considering the 
correlation of the different azimuthal scatters with the total dumpi- 
ness, C, the correlation still exists although being much weaker. It 
is worth to point out also how the trend with the X-ray energy is 
opposite in the two cases. Since clumps embed lower temperature 
ICM, softer bands are more correlated to C with respect to the hard 



Table 1. Spearman's rank correlation (rg) of the values of C and Ctz 
(second and third column, respectively) with the azimuthal scatter in the 
X-ray surface brightness in different bands (second to fourth row, respec- 
tively) and in the y-parameter profiles (fifth row). The values analysed 
are taken from the 30 relaxed haloes considering the bins in the range 
0.15 < r/R 20 o< 1.5. 





c 


Ctz 


ctsb (0.3-0.8 keV) 


0.32 


0.59 


CT SB (0.5-2 keV) 


0.30 


0.60 


<T S B (2-10 keV) 


0.23 


0.61 


(7y 


0.35 


0.68 



one. On the other side, when clumps are excised, the densest gas is 
also hotter thus the correlation Ctz-o~sb increases at higher energy. 

More in the detail, we show in Fig.[8]the correlation between 
Ctz and the scatter in the 2-10 keV band (top panel): although 
there is a significant dispersion, the trend of increasing dumpiness 
with increasing value of <tsb is clear. The same applies also when 
considering the scatter in the y-parameter (bottom panel). When 
analysing the dispersion of the relation, it is clear also the trend of 
higher values of Ctz for increasing radii (indicated by the different 
colors), as already discussed for Figs.|4]and|5] 

Since it is not possible to define a priori the relation between 
scatter and dumpiness, we proceed empirically by introducing the 
following function 



cro ro 



(6) 



that connects the azimuthal scatter a and the distance from the 
center r with an estimate of the dumpiness Ctz"- The choice of 
this expression is made in order to have a simple increasing func- 
tion of both <r and r (which holds for cro > and ro > 0), and 
CfV = 0,r = 0) = l. 

Given the two free parameters, <ro and ro, for every observ- 
able quantity object of our analysis (i.e. azimuthal scatter in the X- 
ray surface brightness and in the y-parameter) we determine their 
best-fit values by fitting the expression of eq. ((6) using the points 
displayed in Fig. [8] i.e. Ctz as a function of a and r. Again, we 
restrict this calculation to the 30 relaxed clusters and to the range 
0.15 < r/i?200< 1.5. The results are shown in Fig. [8] together 
with the residuals, for the 2-10 keV band and for the y-parameter. 
It is clear that the best-fit function provides a good global descrip- 
tion of the relation Ctz-(u, r), with almost all of the points that 
have an estimate within 10 per cent of the true value. Furthermore, 
we note how the diagonal dashed line corresponding to best-fit re- 
lation in the limit of r = defines well the "forbidden" region of 
the Ctz-o plane, below the line itself. Table [2] shows the best-fit 
values for all the three X-ray bands and for the thermal SZ effect. 



5 IMPROVING MASS ESTIMATES 

We have shown in the previous section that it is possible to ob- 
tain an estimate of the large-scale inhomogeneities (Ctz) from the 
azimuthal scatter in the X-ray surface brightness and y-parameter. 



In the fitting procedure we omit to assign errors to our data. 
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Figure 8. Upper panel: correlation between the azimuthal scatter in the 2- 
10 keV band and C'tz. The points considered correspond to the sample of 
the 30 relaxed haloes in the range 0.15 < r/R2oo< 1.5 (a total of 990 
points, 33 for each halo) and their colors indicate the radial distance, in- 
creasing from blue to red. The quantities indicated (<ro and ro in units of 
i?20o) are the best-fit values obtained using the formula of eq. l(6j and the 
lower part shows the residuals with respect to the model in terms of relative 
difference in the dumpiness value. The diagonal dashed line indicates the 
relation Cji= 1 + cr/ao, i.e. the best-fit formula in the limit of r = 0. 
Bottom panel: same as upper panel but for the scatter in the y-parameter. 



The question that arises is whether it is possible to use the observed 
value of a(r) as a function of the distance from the cluster centre, to 
improve the estimates of the gas density profile and, consequently, 
the measurement of Af gas and M-^ c . 



5.1 Correcting the gas density profile 

The estimate of the gas density profile p(r) from the X-ray surface 
brightness is affected by the problem that the corresponding emis- 
sivity ex(r) in a fixed band depends on the squared gas density: 
ex oc < p 2 >v- Given the definition of eq. ((2), an observer that 
ignores the residual dumpiness of the gas and assumes spherical 
symmetry obtains 



and use the relation of eq. ((6) to obtain an estimate of Ctz, we can 
correct our measurement of the density profile as follows 



p(r) 



Px(r) 



(8) 



As examples, we show in Fig.|9]the application of this method 
to the two clusters of Fig. [2] The relaxed D17-a cluster (left panel) 
shows reconstructed density profiles very close to the true valufl 
with differences smaller than 2 per cent up to -R200 ■ On the other 
side, the corresponding uncorrected X-ray profile overestimates the 
true value by 5-10 per cent over the whole r > OAR200 range. 

When looking at the D12-a cluster (right panel), the presence 
of the disturbing structure still produces an overestimate of the cor- 
rected density profiles of about 5-10 per cent in correspondence of 
the dumpiness peak. However, even in this case the improvement 
with respect to the original px profile overestimate is remarkable, 
through all the radial range. 

This method can be directly ap plied to the observed density 
profiles and scatter of lEckert et al.l d2012l) . In fact, they obtained a 
measurement of ass in the soft 0.5-2 keV band in their 31 ROSAT- 
PSPC objects (we refer to their Fig. 7, entire sample) up to i?200- By 
using this quantity to estimate the residual dumpiness, we obtain an 
approximate value of C|f = 1.2 at r ~i?200, that corresponds to an 
overestimate of ~ 8 per cent in the gas density. This consideration 
leads to mitigate the existing conflict between their observed den- 
sity profiles and the simulated ones, although not enough to solve 
the problem completely. 



5.2 Correcting the gas mass bias 

If its gas density profile p(r) is known, then the gas mass of the halo 
enclosed by -R200 can be obtained from the following formula: 



M, 



47T 



-^200 



p(r) r dr . 



(9) 



However, the density bias described by eq. ((7) reflects also 
into the gas mass estimate, 



Mgas.x = 4tt / y/Cn{r)p{r) r 2 dr . 
Jo 



that leads to a bias in the gas mass measurement 



6(M g 



gas,X 



M s . 



(10) 



(11) 



which is always > reaching null value only in the limit of a per- 
fectly uniform medium (Cn= 1). 

Given our set of simulated clusters and groups, we can obtain 
the expected value of 6(A/ gas ) from their density and dumpiness 
profiles. The results are shown in the top panel of Fig. [10] typical 
values correspond to overestimates of 5-10 per cent, with the more 
perturbed haloes that can reach overestimates around 20-30 per 
cent. 

If we substitute in eq. l[9} the corrected density profile of 



px(r) = v^PW , (7) 
to a systematic overestimate. 

If we consider the observed azimuthal scatter profile, cr(r), that small clumps have been efficiently removed in X-ray observations. 



4 Here we are considering the density profiles after applying the volume- 
which corresponds to a systematic overestimate. clipping method described in Section [375] so we are implicitly assuming 
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Figure 9. Sketch of the gas density correction method applied to the two clusters shown in the maps of Fig. [2] D17-a (left panel) and D12-a (right panel). 
We show the radial dependence of the relative differences between the true gas density and the uncorrected measured density (solid line), the density after 
applying the X-ray and SZ corrections (dot-dashed and dotted line, respectively). The horizontal dashed line shows the ideal case of perfect density profile 
reconstruction (Ap = 0). 



eq. l[8j, we obtain a new estimate of the gas mass 

M gas = 4tt / p(r) r 2 dr = 
Jo 

which provides a good approximation of M gas as much as C^* is 
an accurate estimate of the residual dumpiness. 

We verify the accuracy of this method by applying it to our 
simulated haloes. The central and bottom panels of Fig. [10] show 
the distribution of the bias after applying the correction consid- 
ering the scatter in the 2-10 keV band surface brightness and in 
the y-parameter, respectively. We report in Table |2] the values of 
the median of the corrected and uncorrected bias distributions to- 
gether with their dispersion also for the other two X-ray bands, 
for our sample of 30 relaxed haloes. In all cases the systematics 
for relaxed systems is completely eliminated, with the halo to halo 
scatter which is reduced as well. A very small number of highly 
perturbed haloes keep having overestimates of 10-20 per cent: this 
reflects the presence of other types of asymmetries in the profiles 
not directly associated to the dumpiness of the ICM. Nevertheless, 
when calculating the median and the quartiles of the whole sample 
distribution, as shown in Table|3] we can see that our method is still 
effective even if applied "blindly" to all objects. 

In Fig.fTTJwe show the radial dependence of &(Af gas ) and of 
the corresponding bias after applying the X-ray and SZ corrections. 
Even close to the core, the uncorrected bias is significant (~ 5 
per cent) and it grows linearly up to 10 (20) per cent at -R200 for 
relaxed (perturbed) systems. We observe that already at distances 
corresponding to -R500 (~ 0.7i?20o) the expected overestimate is 
already close to the values at -R200, with differences of the order 
of 2 per cent. When applying the corrections, we can clearly see 
that for relaxed objects the systematics is completely removed with 
value a of ± 2 per cent up to 2R200 ■ Even for perturbed haloes 
the improvement is consistent (a factor of ~3) through the whole 
radial range. 



5.3 Correcting the hydrostatic-equilibrium mass bias 

A frequently used method to determine the total mass of a galaxy 
cluster relies on the assumption that the ICM is in hydrostatic equi- 
librium. The measured density and temperature profiles are then 
used to determine the total mass Mh c (< r) enclosed by a given 
radial distance r 

m r) = ^ + dj^rn) t (13) 

G/im p \ alogr alogr J 

where T(r) is the (mass-weighted) temperature profile, and p is 
the mean molecular weight in units of the proton mass, m p . This 
method is known to be affected by several systematics associated 
to the break of the spherical symmetry and to the presence non- 
thermal pressure sources that affect particularly cluster outskirts 
(see Section [TJ. Here, instead, we focus on determining the bias 
associated to the residual dumpiness of the ICM. 

It is easy to show that when accounting for the gas density 
bias of eq. 0, the X-ray measured hydrostatic-equilibrium mass 
Afhe,x formula contains an additional term associated to Cu, thus 
becoming 

M (^,\ fe B T(r)r ^ rilogp(r) dlogT(r) , 

M he ,x(< r) = — — 1 — h 

Gpm p \ alogr alogr 

+ i^M) . (14) 
2 a log r J 

Since Cuij) generally grows with r (see Figs.[4]and[5}, unlike p(r) 
and T(r), this turns into an underestimate of the value of Mh e - To 
study this effect in detail, for every object of our sample we com- 
pute the true value of A/ho and Mh c ,x by adopting the following 
procedure. We first fit T(r) in t he range 0.8-1.2 i?200> with the 
formula of Vikhli nin et al.l (20060 Then, in the same radial range, 
we fit p(r) and the biased density px(r) profiles with a /3-model 

5 Here we consider only the external part of their profile, i. e. the formula 
of their eq. (4), after fixing a = 1 and 6 = 2, thus leaving only three free 
parameters: normalisation, core radius and the external slope c. However, 
since our objective is to determine the bias associated to the uncertainties in 
the density profile, any possible effect introduced by changing the tempera- 
ture fitting procedure is negligible. 
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Table 2. Best-fit values of the parameters ctq and tq (second and third column, respectively, with ro in units of i?20o) of the relation Cp St (<r, r), being <r the 
azimuthal scatter of the surface brightness in the three X-ray bands (second to fourth row) and the scatter of the y-parameter profiles (fifth row). The fourth 
column shows the median (in percent units) of the distributions of bias in the value of M gas at R200 f° r our sample of 30 relaxed haloes, together with the 
difference with respect to the upper and lower quartiles (indicated as error), after the dumpiness correction of eq. (II 21 . The fifth and sixth column show the 
same quantities for the distribution of M^ c and / ga s, respectively. In the last row we report the values corresponding to the uncorrected biases. 
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ro 


6(M gas ) % 


b(M bc ) % 




<tsb (0.3-0.8 keV) 
cr SB (0.5-2 keV) 
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16.02 
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Table 3. Same as Table|2]but for the whole sample of 62 haloes. The best-fit <ro and ro are not quoted here since we assume the same values of the relaxed 
sample. 
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dCavaliere & Fusco-Femianolll978l) and use these results to com- 
pute A'/h c and Mh e ,x. respectively, with eq. d!3t for r =i?2oo- 

The top panel of Fig.[T2]shows the distribution of the expected 
value of 6(Mhc) for our sample of simulated clusters and groups. 
While it is confirmed that for the majority of the cases (45 over 
62) the residual dumpiness introduces an underestimate, since the 
clumping factor of single objects can be affected by local variations 
(i.e. Ciz(r) is not perfectly monotonic), we end up with about a 
quarter of our sample in which the measured Mh c is higher than 
the true value. As a general remark, the median effect of about —2 
per cent is small when compared to the other systematics, and with 
respect to the intrinsic dispersion (~ 6 per cent difference between 
the upper and lower quartile). We note, however, that for the most 
disturbed systems the error can reach 10 per cent or more in either 
direction. 

We repeat the same procedure done for 6(M gas ) (see Sec- 
tion \52& to see whether our clumpiness-correction method is ef- 
fective in reducing the value of b(M bc ). The results for the three 
X-ray bands and for the SZ effect are shown in Table [2] and [3] and 
in the central and bottom panel of Fig.[l2]for the hard 2-10 keV 
band and for the SZ effect only, respectively. Also in this case we 
observe a global improvement of the measurements, with the sys- 
tematics that disappear in all cases. On the other side, the halo to 
halo scatter is not significantly reduced. 



5.4 Correcting the gas fraction bias 

The two biases associated to the residual dumpiness previously de- 
scribed (positive bias for M gas and negative for A/h c ) add together 



when measuring the gas fraction / gas . In fact, being the measured 
value 
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/gas.X 



ga s,X 



Mi 



(15) 



hc,X 



the corresponding bias is 



s _ /gas,X ~ / gaB ft(Af gM ) ~ b(M hc ) 

b(/sas) = / gas = i + w • (16) 



When computing the value of i>(/ gas ) at r —R200 for our sam- 
ple of 62 haloes, we can see from the distribution shown in the 
top panel of Fig. [T3] that this results in a systematic overestimate 
on average of about 10 per cent. The ICM inhomogeneities cause 
also a large spread in the measured value, with 1 1 objects having 
6(/ gas ) > 0.2 (three are outside the plot range). 

We verify the efficiency of our bias-correction method for / gas 
by computing the expected value of b(/ gas ) once the corrected val- 
ues of 6(M g as) and b(Mh c ) (see Section [5T2l and [531 respectively) 
are used in eq. l !16t . We show the results in Table [2] and [3] and in 
the central and bottom panel of Fig. |T3] (2—10 keV band and SZ 
effect only, respectively). Again, in all cases the average bias is sig- 
nificantly reduced, with <tsb f° r the 2-10 keV band and a y that 
provide the best results. In this case, as for b(M^ c ), the applica- 
tion of our method allows also to reduce the scatter in the measured 
values. 
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Figure 10. Top panel: histogram of the estimated gas mass bias associated 
to the residual dumpiness for our sample of 62 clusters. Perturbed haloes 
are marked by the shaded areas. Central and bottom panels: same as top 
panel but after the correction obtained by estimating the value of C'-jz from 
the scatter in the 2-10 keV surface brightness and y-parameter, respectively. 



6 SUMMARY AND CONCLUSIONS 

In this paper we have analyzed a set of 62 simulated clusters and 
groups, obtained with different physical prescriptions, with the ob- 
jective of providing a global characterization of their density in- 
homogeneities in the regions close to the virial radius. We have 
described a method that allows to separate the ICM dumpiness as- 




Figure 11. Bias in the value of M gas as a function of distance from the 
centre for our reference physical model without correction (solid lines) and 
with the correction obtained by estimating the value of C-jz from the scatter 
in the 2-10 keV surface brightness and y-parameter (dotted and dot-dashed 
lines, respectively). The black lines represent the median of the 30 relaxed 
objects and the grey-shaded regions enclose the quartiles of the uncorrected 
and of the X-ray corrected biases (the quartile region of the SZ corrected 
bias has a similar size). The red lines represent the median value for the 
perturbed sample. The horizontal dashed line corresponds to the value of 
6(M gas )=0. 



sociated to small clumps from a "residual" one that corresponds to 
large-scale inhomogeneities, and discussed how the latter depends 
on the mass and the dynamical state of the halo. Finally, we have 
discussed how the presence of large-scale inhomogeneities can bias 
the estimates of M gas , Mh c and / gas and provided a method to re- 
duce this bias by using a directly observable quantity: the azimuthal 
scatter in the X-ray surface brightness profiles or in the thermal SZ 
ones (y-parameter profiles). 

Our results can be summarized as follows. 

( i) As expected, the degree of global dumpiness in our simulated 
objects depends mainly on the presence/absence of radiative cool- 
ing, making about one order of magnitude difference. Once cooling 
is included, additional feedback mechanisms do not change signif- 
icantly the clumping factor. 

(ii) When considering only the contribution of emitting gas, we 
obtain values of Cc± 2 — 3 at -R200 ■ When compared to the value of 
~ 16 claimed bv lSimionescu et all J201 lh . we conclude that their 
estimate can not be reproduced by our models, even neglecting the 
possibility of identifying emitting clumpy structures. 

(Hi) We introduce the concept of residual dumpiness, Cn, to 
quantify the amount of large-scale inhomogeneities (departure 
from spherical symmetry, presence of filaments), that corresponds 
to the bulk dumpiness once obvious bright condensed clumps are 
masked out. This quantity is independent of the physical model 
assumed, while it is sensitive to the dynamical state of the halo: re- 
laxed objects have Cn~ 1.2 at -R200, while dynamically perturbed 
ones have on average C'tz^ 1.5. 

(iv) The residual dumpiness of the ICM causes a significant over- 
estimate in the measurement of M gas from X-ray observations of 
the order of ~ 5 — 10 per cent for the more relaxed objects up to 
~ 20 — 30 per cent for the more perturbed ones. 

(v) A smaller negative bias of about 2 per cent is present also 
in the measurement of Mho- Consequently, the combination of the 
two biases results in an overestimate of / gas , with average values 
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Figure 12. Same as Fig. llOl but for the bias in the hydrostatic-equilibrium 
mass. One perturbed haloes, D03-fl, is not shown since its values are outside 
the plot range being, from top to bottom, b(iWj, e ) = —0.22, —0.21 and 
—0.19, respectively. D27-a is also outside the plot range in the top and 
central panels with &(M[ 1C ) = —0.25 and —0.17, respectively. 
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Figure 13. Same as Fig. llOl but for the bias in the gas fraction. Three per- 
turbed haloes, D03-a, D18-C/ and D27-0, are not shown in the top panel 
since their values are outside the plot range: 6(/ gas ) = 0.35, 0.33 and 
0.65, respectively. D27-a is outside the plot range also in the central panel, 
having f>(/ gas ) = 0.30. 



of ~ 10 per cent, and an intrinsic scatter of ~ 9 per cent (interquar- 
tile range). These biases are lower when compared to other known 
systematics. 

(vi) The residual dumpiness correlates well (rs = 0.6 — 0.7) with 
the azimuthal scatter of the X-ray surface brightness and of the y- 
parameter profiles. This allows us to obtain an analytical formula to 



estimate it as a function of two observable variables: the azimuthal 
scatter and the radial distance. 

(vii) This relation provides a method to correct the gas density es- 
timates, making it possible to improve consistently the accuracy 
of the Afg as measurements. With this method the systematics de- 
scribed above disappears completely for relaxed haloes from out- 
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side the cluster core up 2i?2oo- For perturbed clusters/groups the 
overestimate is reduced by a factor of about 3. 
f viii) Finally, this method works also in eliminating the bias associ- 
ated to the measurement of M hc and / gas . However a large intrinsic 
scatter (5-7 per cent, in terms of interquartile range) between the 
different objects remains. 

Overall, our results show how the study of the outskirts of 
galaxy clusters and groups is important for the measurement of the 
gas mass and gas fraction, and how the combination of simulations 
and observations can improve their precision. A possible extension 
and improvement of this analysis may be investigating the correla- 
tion of Cn with other inhomogeneities parameters such as the halo 
ellipticity, or by determining how the Cn(p, r) relations may vary 
as a function of the observed relaxation parameters of the haloes. 
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APPENDIX A: INDENTIFYING CLUMPS WITH 
OBSERVATIONS 

Our volum e-clipping metho d, described in Section |3~l2l and more 
in detail in lRoncarelli et alj d2006bl) . allows us to separate between 
the gas belonging to clumps and to the bulk of the cluster based on 
theoretical considerations. Here we provide more information on 
the physical properties of these clumps and investigate up to which 
extent they may be detected in simulated X-ray surface brightness 



10 1D 




1 2 1 i i i i i i 

10° 10' 10 2 1 3 10* 1 5 10 6 1 7 

P/Pb 

Figure Al. Phase diagram of the SPH particles outside O.5.R200 of the D17- 
a simulated cluster for our reference model. The black points indicate the 
particles of the bulk. The particles identified as clumps by our algorithm 
are marked in green if they fall inside a contour region of the right panel of 
Fig. lA2l and red otherwise. 

maps, and whether the different methods would introduce signifi- 
cant biases in the Cn determination. 

We show in Fig. lAll the T — p scatter plot of the SPH parti- 
cles in the outskirts of the D17-a relaxed cluster. The amount of 
gas that our algorithm associates to the clumpy phase (red and 
green points) corresponds to about the 3 per cent of the total gas 
mass inside R200 for this halo. We verified that in perturbed sys- 
tems it is slightly higher (4—5 per cent). By analysing the plot, the 
multi-phase nature of these clumps shows up clearly. Part of the 
gas is associated to the cold star-forming phase at p > 10 3 pb and 
T — 10 4 — 10 5 K. These particles are responsible for most of the 
global ICM dumpiness (see the comments on Fig. [4]) but would not 
influence any X-ray measurement since their temperature is too low 
to produce any significant emission. The majority of clump parti- 
cles, which usually embed the previous ones, are instead at higher 
temperature (T > 10 6 K) and can eventually produce a detectable 
X-ray emission. 

The maps of Fig. I A2 1 show the 0.5-2 keV surface brightness, 
up to 2i?200 of the bulk (left panel) and of the clumps (right) of 
the D17-a relaxed cluster, as defined by our method. The great ho- 
mogeneity of the bulk map shows clearly that no clumpy structure 
is missed by our filtering method and that any possible detectable 
inhomogeneity must necessarily be associated to the clump phase. 
The possibility of detecting them depends on how much their signal 
is brighter with respect to the bulk one. To this purpose we identi- 
fied the map regions where the signal-to-noise ratio exceeds a value 
of 3, by considering the bulk map as a reference for the noise. Most 
of the clumpy structures are at least partially identified with this 
method, with an increasing efficiency towards the outskirts. 

In Fig. IA3I we show the surface brightness of the identified 
clumps as a function of distance from the centre, compared to the 
average cluster surface brightness profile. Since the central part of 
these clumps is usually very bright, in most of the cases their aver- 
age surface brightness exceeds the cluster one by more than 5. It is 
clear also, that these structures can be identified also outside -R200, 
since their signal is well above the unresolved X-ray background, 
at least in the case of a massive halo. 

As a drawback of our method, when in a given radial bin no 
clear clumpy structure is present we consider as clumps also a small 
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Figure A2. Maps of the soft (0.5-2 keV) X-ray surface brightness of the D17-a cluster for our reference model in units of counts s — 1 cm -2 arcmin - 2 . Each 
map is centered on the cluster centre and encloses a square of 4i?200 per side. Left panel: the gas which is identified as 'bulk' by our theoretical filtering 
method (99 per cent volume). Right panel: the gas which is identified as 'clumps' (one per cent volume). The contours enclose the regions where the projected 
clumps signal is more than 3 times the bulk one, mimicking a S/N>3 threshold. 



fraction of diffuse gas whose signal is not high enough to be iden- 
tified. This happens in particular in the regions close to the centre 
where the ICM is more uniform, as it can be seen both in the right 
panel of Fig. lA2l and by the red points of Fig. IA2I A precise de- 
termination of the fraction of clumps that would be missed by real 
X-ray observations clearly depends on instrumental details and is 
beyond the aim of our work which focuses on the large-scale in- 
homogeneities. However, using our mock maps we can provide an 
estimate of the amount of missed clumpy gas, together with its in- 
fluence on the bias of the gas mass measurements. We verified that 
the detected clumpy ICM (e.g. green particles in Fig. IA2t is the 
~20 per cent of the total amount of gas associated to clumps by 
our volume-clipping method. In order to have an estimate of its 
impact on the mass gas measurement, we repeated our analysis by 
reducing by a factor of 5 the volume threshold (i.e. the 0.2 per cent 
of the volume of each bin). We obtained a value of fo(M gas ) = 7.7 
per cent for the relaxed sample, with respect to the previous 6. 1 (see 
Table[2} obtained with our more conservative threshold. Therefore, 
given the relatively small difference, we conclude that the impact 
on our results of any unresolved clumpy gas is minor with respect 
to the one associated to large-scale inhomogeneities. 
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Figure A3. Soft (0.5-2 keV) X-ray surface brightness of the clumps (green 
diamonds) identified with the method shown in Fig. IA2I as a function of 
distance from the centre. The black solid line represents the average surface 
brightness of the cluster bulk. The red dashed lines indicate the measure- 
ments of the unresolved X-ray background in the same band from the COS- 
MOS survey jElvis et alj 2009. ~200ks exposure, uppe r line) and from the 
4Ms Chandra Deep Field South ICappelluti et alj2012l lower line). 



APPENDIX B: RELAXED AND PERTURBED HALOES 
FROM AN OBSERVATIONAL POINT OF VIEW 

In Section 13.21 we have described our classification of simulated 
haloes into relaxed and perturbed according to the robustness of the 
determination of Ciz. Here we explain how this criterium, which is 



based on purely theoretical considerations, actually matches well 
with a possible classification done with observational methods. 

To this purpose we have created a 0.5-2 keV band sur- 
fac e brightness map for each of our simulated objects. Follow- 
ing [CassantTetal] d2010h . we have then computed, over the radial 
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Figure Bl. Power-ratio P3/P0 versus the centroid shift w estimated from 
the 0.5-2 keV surface brightness maps in the radial range 0.1-1 R200 f° r 
the sample of our 62 simulated clusters and groups. The two objects indi- 
cated with arrows have values outside the plot range. The two dashed lines 
indicate the limits of w = 0.025 and P3/P0 = 10~ 7 that provide the best 
match between the two definitions of relaxed/perturbed halos (see text). 



range 0.1 — 1-R2oo> the centroid shift w, defined as the standard 
deviation, in units of -R200, of the projected separation between 
the X-ray peak an d the centroid, and the power-ratio P3/P0 (see 
lBuote&Tsailll995h . that is the lowest normalized moment of the 
X-r ay surface brightness c learly connected to substructures (see, 
e.g jBohringer et al J2010h . We show in Fig. IB li the position of our 
simulated systems in the P3/P0 vs w plane. Haloes defined as re- 
laxed according to our definition (marked in black) have on average 
lower values of all the parameters, thus sitting on the lower left cor- 
ner of the plot, while perturbed ones tend to show a high value in 
at least one of the two parameters. 

The dashed lines (10 = 0.025 and P 3 /P = 10~ 7 ) show the 
limits that roughly correspond to the best match between the two 
classifications. When defining observationally the relaxed haloes 
as the ones laying on the bottom-left quadrant of the plot, and per- 
turbed otherwise, we end up with 53 over 62 haloes matching the 
corresponding theoretical definition. 



